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Abstract - We present numerical evidence that the dynamical zeta func- 
tion C,{z) and the Fredholm determinant F(z) of intermittent maps with a 
neutral fix point have branch point singularities at z = 1. We consider the 
\Q \ power series expansion of l/£(z) and F{z) around z = with the fix point 

pruned. This power series is computed up to order 20, requiring ~ 10 5 peri- 
odic orbits. We also discuss the relation between correlation decay and the 
nature of the branch point. We conclude by demonstrating how zeros of zeta 
functions with thermodynamic weights that are close to the branch point can 
be efficiently computed by a resummed cycle expansion. The idea is quite 



similar to that of Padee approximants, but the ansatz is a generalized series 



expansion around the branch point instead of a rational function. 



In this note we will present numerical evidence that dynamical zeta functions and 
Fredholm determinants exhibit branch point singularities for intermittent maps. We will 
consider a family of one dimensional maps x i— > f(x) on the unit interval, with 

( ,_jx + 2 s x 1+s 0<x<l/2 
na;j ~| 2a? -1 1/2<x<1 ' [) 

where s > 0. For s = the map would just be the binary shift map, which is uniformly 
hyperbolic, but for s > it is intermittent; the fix point x = is neutrally stable: 
/'(0) = 1. The map admits a binary coding, we associate the letter with the left leg, 
and 1 with the right leg. The neutral fix point now corresponds to the periodic orbit 0. 
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The dynamical zeta function is denned by 

i/CW=Il( 1 - w\) • ( 2 ) 

p \ p\ 

The product in (|2|) runs over all primitive periodic orbits p, having period n p and stability 
A p = ^-\ x=Xp with x p being any point along p. 

In our considerations we will prune the neutral fix point 

1/C(*) = (1 - z)l/l{z) , (3) 

where 

i/CW = Il( 1 -nj) • W 

p^o 1 pl 

We will also consider the Fredholm determinant 

n*) = n n a - iaia^) • (5) 

In the previous case, it is not essential to prune the neutral periodic orbit, but when 
considering the Fredholm determinant it is essential, the extra factor (1 — z)°° would of 
course be devastating for our investigations. 

We will consider the series expansion of the functions \/Q{z) and F(z) 

oo 

V(» = l-E^ > (6) 
1=1 

oo 

F(z) = l-£M* . (7) 
i=l 

The nature of the leading singularity will be reflected in the asymptotic behaviour 
of the coefficients of these power series. To get an idea what this asymptotic behaviour 
may be we consider the fundamental part of a cycle expansion of 1/C(z) |j 

oo n+1 

VC»«i-£t= • (8) 

n=0 A 0"l 

The dependence of the stabilities A^^- on n can be estimated by replacing the difference 
equation x n+ i = x n + 2 s x^ +s by the differential equation^] 



(in 

having solution 



2 S < +S , (9) 
xn = [xq s ~ s2 s n}- l ' s . (10) 
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First, because x n ~ 1 we get a relation between xq and n 

xo-n- 1 /* , (11) 

telling us how deep into the intermittent region the orbit n l penetrate. We are only 
interested in the dependence on n and have discarded other information. We can now 
obtain the stabilities by differentiation 

Acrn = 2^ - 2x s n + % (S+1) ~ n^' . (12) 

This suggest that Z(z) contains a singularity of the type 

{l-z) l/s l/s 
(l-z) 1 /«l g(l-z) l/seN+ , { } 

as can be realized through the Tauberian theorems for power series. 

A somewhat more refined, although closely related estimate is obtained by the piece- 
wise affine approximation of ref. ||, leading to the same predictions. In j|] there are 
also much stronger arguments that z = 1 is a branch point but it is still not proven that 
l/£(z) can be analytically continued outside the unit disk. That this can actually be 
done is suggested by the numerical data we are now going do discuss. 

In order to numerically determine the coefficients {aj} we simply expand the product 
(0) using all periodic orbits up to period 20. This set contains 111011 periodic orbits 
which means that we cannot go very much further within reasonable computer time. 
When expanding the Fredholm determinant (|J) we use Euler's formula to expand the 
inner product: 

m = n n a - wx^) ( i4 ) 

j(j-i) 
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K U |A P P'lTLi(l-A- fe ) 

Periodic orbits are determined by a Newton-Raphson procedure. To this end we look for 
fix points of some iterate of the inverted map, choosing branch according to the symbol 
code. 

The coefficients for three parameter values are plotted in figs. 1-3 together with the 



expected slope according to eqs. (|8|) and (12|). Our set of data is consistent with the 



expected leading singularity (|T^) and indicate no other singularity close to the unit circle; 
this would have induced oscillations superposed on the sequence of coefficients. 

For the binary shift map (s = 0) we would have a% = l/2\ For a slightly higher value 
of the parameter s = 0.1 (fig. 1) the coefficients for small i conform with this behaviour 
but eventually they bend off, approaching the expected slope. 

One may also observe that whereas eq. (||) seems to provide the correct order of the 
singularity it does not predict the correct size (prefactor) of it. The curvature corrections 
of ref. p| exhibit the same type of singularity as does the fundamental part. 
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Why is this interesting? It is known that the spectra of zeta functions and Fredholm 
determinants are, at least in some cases, related to the (typical) decay of correlations. 
For exponentially mixing system the mixing rate is given by the size of the gap between 
the leading and next-to-leading zeros of F(z) We believe that a similar coupling can 
be made system without a gap. Let us take the (formal) trace of the transfer operator 



1 oo 



5n,rn p 1 f -n F '( z 



trC n = 5(x- f n (x))dx = V n p V r—^-~: = — / Z^-^-dz (15) 



P ^\det(l-A;)\ 2iriJ lzl<1 F(z 

where e is a small number. The trace is formal because it makes no explicit reference to 
eigenvalues of an operator, it is just the trace over the integral kernel of the operator. 
We claim that this trace serves as an archetype correlation function. If there is a gap, 
residue calculus tells us that the trace will approach unity exponentially fast and the rate 
is provided by the size of the gap. But could this really work for the intermittent case 
where the leading zero is connected by a branch cut running along the line Im(z) = 
Re(z) > 1. Let us assume that F(z) is holomorphic and zero-free, except along the cut, 
in the disk 1 < \z\ < C where C > 1. The value of the (15) for large n above will 



be governed by the vicinity of z = 1 and can be evaluated asymptotically. Inserting F 
with the suggested branch point singularity into (15) we get the following asymptotic 
behaviour for the trace 

1 + C/n 1 / 3 - 1 < s < 1 
trC n ~ { 1 + C/logn s = l (16) 

1/S 8<1 

For < s < 1 this suggests that typical correlations should decay as ~ l/n 1 / 8-1 which 
indeed agrees with the rigorous results || . The failure of the trace to approach unity for 
s > 1 reflect the fact that the invariant density is not normalizable anymore. 

There are indications that the identification between the behaviour of the formal 
trace and the typical correlation functions is possible also for the Sinai billiard which 
seems to exhibit the decay law C(t) ~ 1/t B |. 

Our findings are also interesting from a more practical point of view. Consider the 
problem of computing the topological pressure which amounts to compute the leading 
zero of 

i/c>,/?) = ri( 1 -^) • 

If is close to, but less than unity this leading zero will be close to z = 1 and it will 
be extremely inefficient to compute it from the truncated power series in z. It is natural 
to try to expand l/£(z,/3) in a generalized power series around z = 1. If the leading 
singularity is of the form (1 — z) a the simplest possible expansion would be 

oo oo 
1/C(Z, 0) = - zY + (1 - z) a zf . (18) 

i=0 i=0 
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According to our previous findings we expect that a = (3(s + l)/s — 1. Suppose now 
that we replace these infinite sums by finite sums of increasing degrees, n c and n<2, an d 
require that 

$>(1 - zf + (1 - zY £ *(1 - *Y = I/O, P) + 0(z n+1 ) . (19) 

i=0 j=0 

If ric + n^ = n + 1 we just get a linear system of equation to solve in order to to determine 
the coefficients q and di from those of the series expansion around z = 0. This strategy is 
quite similar to that of Pade approximants. It also natural to require that \rid+oi— n c \ < 1. 
So far we have assumed that a is an non integer. The case with integer a can be worked 
out in close analogy. 

To test the idea we choose (arbitrarily) the parameters s = 0.7 and f3 = 0.9. In fig 
4 we plot the leading zero versus truncation length n determined from expansion ( |i~9|) 
and the expansion around z = 0. The improvement is obvious. However, we do not 
claim that the simple expansion (|l~8l ) is entirely able to capture the complicate analytic 
structure around z = 1. 

I would like to thank Viviane Baladi for gently persuading me to undertake this study. 
This work was supported by the Swedish Natural Science Research Council (NFR) under 
contract no. F-FU 06420-303. 
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Figure 1: Expansion coefficients for the functions Z{z) and F(z) for the parameter 
values s = 0.1 
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Figure 2: Same as fig 1 but with s — 1. 
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Figure 3: Same as fig 1 but with s = 2. 
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Figure 4: Leading zero for s = 0.7 and (3 = 0.9 versus truncation length determined 
from the generalized series expansion around z — 1 and the power series around 
z = 
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